Origin of the "waterfall" effect in phonon dispersion of relaxor perovskites 
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We have undertaken an inelastic neutron scattering study of the perovskite relaxor ferroelectric 
Pb(Zn 1/3 Nb 2/ 3)03 with 8% PbTi0 3 (PZN-8%PT) in order to elucidate the origin of the previously 
reported unusual kink on the low frequency transverse phonon dispersion curve (known as "waterfall" 
effect). We show that its position (qw) depends on the choice of the Brillouin zone and that the 
relation of qw to the size of the polar nanoregions is highly improbable. The waterfall phenomenon 
is explained in the framework of a simple model of coupled damped harmonic oscillators representing 
the acoustic and optic phonon branches. 
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I. INTRODUCTION 



Ferroelectric perovskites Pb(Zn 1 / 3 Nb 2 /3)03 (PZN) 
and Pb(Mg 1/3 Nb2/3)03(PMN) and related materials 
have recently attracted a great attention, among oth- 
ers, due to the excellent piezoelectric properties 1 of their 
solid solutions with PbTiC>3. Their average ABO3 per- 
ovskite structure is perturbed by a short-range occupa- 
tional ordering on B-sitc positions!^ It is known that 
this nanoscopic inhomogcneity is responsible for relaxor 
properties of these materials, such as smearing and fre- 
quency dependence of the dielectric anomaly associated 
with the ferroelectric ordering 4 :, and for other phenom- 
ena reflecting the presence of polar nanoregions (PNR) 
persisting hundreds of Kelvins above the temperature of 
dielectric permitivity maximum*. Unfortunately, it is 
not easy to deal with these inhomogeneities in a sim- 
ple enough model, and the understanding of the relaxor 
properties is not yet satisfactory. 

Since these materials are available as large single crys- 
tals, inelastic neutron scattering can be used to study the 
phonon dispersion curves. The spectrum of the lowest 
frequency transverse excitations usually consists of in- 
teracting and mutually repelling transverse acoustic and 
transverse optic phonon branches. It is known that upon 
approaching a displacive ferroelectric phase transition, 
the lowest zone center optic mode (soft mode) tends to 
zero frequency. Instead of that, it was repeatedly found 
that near the " smeared" phase transition of relaxor ferro- 
electrics, the upper branch appears to drop precipitously 
into the lower branch at a finite value of momentum 
transfer of the order of q w [ — 0.2 A^ 1 . For explanation of 
this "waterfall" effect, Gehring, Park and Shiranc 6 pro- 
posed that the characteristic wave vector q w { corresponds 
to the size of the PNR's, and that due to these PNR's the 
soft branch phonons with q < q w { cannot effectively prop- 
agate in the crystal. The idea was subsequently specified 
within a mode-coupling models assuming a sharp in- 



crease of transverse optic branch damping (inverse life- 
time) for q < q w {. 
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FIG. 1: Low frequency phonon modes in PZN-8%PT mea- 
sured by constant energy scans at 500 K. Measurements were 
performed in (a) 020 Brillouin zone (Q = (q20) scans) and in 
(b) 030 Brillouin zone (Q — (g30) scans). 

Within the past three years, it was established 
that the "waterfall" effect is common to a number 
of relaxor ferroelectrics, including PZN££, PMN 1 ?! 11 ' 12 , 
Pb(Zn 1/3 Nb 2/3 )03 with 8% PbTi0 3 (PZN-8%PT)&, 
PZN-15%PTi£, and PMN-20%PTi4 However, we have 
noticed that the sharp drop (near q rj 0.2 A -1 ) of the dis- 
persion curves of the single-domain tetragonal BaTi03 
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drawn in Fig. 6 of Ref. also strikingly closely reminds 
those of relaxors (compare, for example, with Fig. 2 of 
Ref. 6) . Obviously, the explanation based on the charac- 
teristic size of PNR's cannot be valid for an experiment 
performed on a classic single-domain ferroelectric crystal. 
This letter brings new inelastic neutron scattering data 
showing that such explanation does not hold for relaxors 
either. We demonstrate that the "waterfall" effect can 
be ascribed to an entirely classic interference of lineshape 
anomalies due to the coupled acoustic and optic branches 
without any ad-hoc assumption of an anomalous increase 
of the bare optic branch damping below q w {. 

For the experiment we have chosen a 3.7 g single crys- 
tal of PZN-8%PT, used already in our previous neutron 
study 1 —. This as- grown, optically transparent, yellowish 
single crystal was produced by the high temperature flux 
technique at Materials Research Institute, Pennsylvania 
State University. This promising PZN-8%PT material, 
showing both relaxor and giant piezoelectric properties, 
has been studied recently by a number of techniques, and 
it is the one for which the possible relation of q w f to the 
size of PNR's was firstly invoked^ 

The present experiment was carried out on the re- 
cently upgraded IN8 thermal neutron spectrometer at 
the ILL high flux reactor. The instrument was oper- 
ated with a fixed wave-number of scattered neutrons 
fcf = 3 A -1 , using horizontally focussing (vertically flat) 
Si crystals as monochromator and analyzer. The sam- 
ple was wrapped in a thin Nb foil and mounted in a 
vacuum furnace allowing to reach temperatures up to 
1200 K with a stability better than IK. The sample was 
oriented with the cubic [001] axis vertical. This geometry 
allowed us to explore the (hkQ) scattering plane. Both 
transverse and longitudinal momentum resolution widths 
measured on the (030) Bragg reflection were better than 
0.07 and 0.06 of reciprocal lattice unit (rlu, where lrlu is 
lixja w 1.55 A -1 ), respectively. This is about equivalent 
to resolution widths in a setup with flat PG (pyrolitic 
graphite) crystals and 40' Soller collimators used in the 
previous studies by other authors. The present setup 
using elastically bent Si crystals^ permits to avoid spu- 
rious effects from tails of reflection curves inherent to PG 
mosaic crystals and provides an improved energy resolu- 
tion of 0.8 meV (FWHM), as given by the energy profile 
of elastic incoherent scattering from our crystal and con- 
firmed by an independent measurement on a vanadium 
reference sample. 

The typical constant energy transverse scans in the 
020 and 030 Brillouin zones at 500 K are displayed in 
Fig. 1, and the positions of their maxima arc plotted as 
solid symbols in Fig. 2. The apparent dispersion curve 
obtained in this way from the 020 Brillouin zone data has 
an almost vertical rise up near q w { = 0.13 rlu. This value 
corresponds exactly to the q w { obtained from the 220 zone 
data at the same temperature and same crystallographic 
direction for this material in Ref. |(J 

On the other hand, the position of the upward rise 
observed by us in the 030 Brillouin zone is at a signifi- 
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FIG. 2: Position of maxima in constant energy transverse 
scans performed in 2 (solid squares) and 3 (solid di- 
amonds) Brillouin zones (as in Fig. 1). Maxima are forming 
two distinct "dispersion curves". The inflection points near 
9wf = 0.13 rlu and 0.06 rlu, respectively, are discussed in the 
text. Open symbols show maximum intensity positions in 
constant-g scans. Lines are only guides for eyes. 

cantly smaller value of q w t w 0.06 rlu. As the same [010]- 
polarized modes with wave vectors along [100] are probed 
in both Brillouin zones, the value of q w { can hardly pro- 
vide a measure of the size of PNR's&. On the contrary, 
the remarkable dependence of q w { on the chosen Bril- 
louin zone proves that the observed dispersion is merely 
an apparent dispersion. In fact, it is known that maxi- 
mum intensity positions for response function of coupled 
damped harmonic phonons can be sensitive to the dy- 
namical structure factors (chosen Brillouin zone)^ Fur- 
thermore, the striking difference between the 030 zone 
value of q w { and almost the same values of q w f in 020 and 
220 zones indicates that the waterfall effect is related to 
dynamical structure factor of bare acoustic branch, which 
is known to be very small in the 030 zoneiii*i£ 

As a matter of fact, it turns out that the "waterfall ef- 
fect" can be reproduced in a quite simple model. Let 
us consider the standard model22*2ii22i of two coupled 
damped harmonic oscillators, defined by a 2 x 2 dynam- 
ical and damping matrices D g and T q 

D « " v A(g)* c4 (g) ) ' ^ 

r _ ( r TA (g) T AO (q) \ (r) , 
^ - [r AO (q) T TO (q) ) ' (A > 

where lota(q) and wto(<7) describes dispersion of bare 
acoustic and optic branches, A(g) describes their mutual 
bilinear interaction, Tta^) and Fto(<?) stands for bare 
mode frequency independent damping, Tao (<?) stands for 
viscous interaction (bilinear in time derivatives of bare 
mode coordinates) and q is the reduced phonon wave 
vector in a chosen direction (in this letter, along the cu- 
bic [100]). In the high temperature limit (hut << kT), 
the corresponding inelastic neutron scattering intensity 
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is proportional toi&^i 

J(w, q) = fcT c- 1 f (<?)* . Im [G(u;, <?)] . f(q) , (3) 

G(w,g) = (D,-i.a;r g -w 2 E)-\ (4) 

where E is a 2 x 2 unit matrix and f(q) is a 2-componcnt 
vector composed of dynamical structure factors of bare 
acoustic and optic branches, respectively. At small fre- 
quencies, it is expected that w.|Fao(<?) << so that, 
as e.g. in the Ref.llSl Tao('z) = will assumed in the fol- 
lowing. The wave vector dependence of dynamical matrix 
D g can be conveniently cast in the form which appears 
in the simplest nearest-neighbor interaction models with 
two degrees of freedom per unit cell 2 ^* 2 ^ 2 ^: 

t4 A (g) = Aam 2 (wq/2), (5) 
<4o(«) = c + Bsin 2 (7rg/2), (6) 
A q = dsin 2 (7rq/2), (7) 

where the real parameters A,B,c and d define the 
strength of dispersion of bare acoustic and optic 
branches, the square of soft mode frequency and the 
strength of mode mixing at zone boundary, respectively. 
Let us note that the dynamical matrix parametrized in 
this way shows the same asymptotic behavior for q — > 
as the standard model of Ref. 1261 and at the same time 
provides the correct dispersion behavior also at the zone 
boundary. In the same spirit, we introduce 

r TA (<z) = gsinVq/2), (8) 

r TO (<z) = h, (9) 

where the real parameters g and h define the damping 
of the bare acoustic branch (with required asymptotic 
q 2 dependence) and the q- independent damping of bare 
optic mode, respectively. 

Numerical calculations were performed for a set of real- 
istic values (A = lOOmeV 2 , B = 150meV 2 , c = 15meV 2 , 
d = 100 me V 2 , g — 6meV, h = 6meV). Constant energy 
profiles calculated for the case of equal bare acoustic and 
optic structure factors (f = (0.5, 0.5)) are shown in Fig. 3, 
while phonon dispersion curves obtained by diagonaliza- 
tion of the dynamical matrix D g are shown in Fig. 4. 
Fig. 4 also shows that the positions of the maximum in- 
tensity in constant energy profiles shown in Fig. 3 forms 
an apparent dispersion branch which follows the acoustic 
branch for small q, but then it sharply rises towards the 
upper phonon branch near 0.1 rlu. The shape of this ap- 
parent dispersion curve indeed corresponds well with the 
experimental findings shown for example in Fig. 2. 

As expected, the "waterfall" region of the apparent 
dispersion curve depends noticeably on the dynamical 
structure factor. For example, a decrease of the weight 
of the bare acoustic structure factor shifts the "water- 
fall" position towards the zone center, as it follows from 
Fig. 4. A similar tendency is found in the experimental 
data taken in the "predominantly optic" 030 Brillouin 
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FIG. 3: Constant energy scans at 1-9 meV simulated from 
the model described in the text, for structure factors f = 
(0.5,0.5). 
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FIG. 4: Dispersion curves calculated from the model de- 
scribed in the text. Full and dashed lines correspond to fre- 
quencies of repelled and bare phonon branches obtained as 
square roots of eigenvalues and of diagonal elements of the 
dynamical matrix D 9 , respectively. Full squares (on the right 
side) stands for positions of maximum intensity in constant 
energy profiles shown in Fig. 3, calculated for balanced struc- 
ture factors f = (0.5,0.5). Full diamonds (left side of dia- 
gram) are calculated in a same way for the case of a smaller 
bare acoustic structure factor (f = (0.25, 0.75)). 



zone (Fig 2). In the limit of zero bare acoustic structure 
factor (f = (0, 1)), the waterfall phenomenon completely 
disappears (Fig. 5). On the other hand, with a slight 
change of parameters, the apparent waterfall dispersion 
curve can attain even a negative slope. These results 
show that constant energy scan technique for investiga- 
tion of coupled acoustic-optic branches should be used 
with certain caution. More detailed discussion of the 
model properties and quantitative analysis of the exper- 
imental data is beyond the scope of this letter and will 
be published elsewhere. 

In conclusion, we argue that the "waterfall" kink seen 
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FIG. 5: Dispersion curves calculated from the model de- 
scribed in the text. Symbols have the same meaning as in Fig 
4.; full squares are calculated from constant energy profiles 
calculated for the case of zero bare acoustic structure factor 
(f=(0,l)). 



constant energy sections of scattered intensity, while no 
such kink is present in the dispersion of bare or coupled 
branches themselves. In the light of these results, there is 
no reason for the position g w f of the waterfall kink to be 
related to the size of the polar nanoregions. First of all, 
there is no need to assume the unusual increase of bare 
optic branch damping below the waterfall wave vector g w f 
in the present model. Secondly, the "waterfall" position 
is very sensitive to the dynamical structure factor both 
in our model and in our experimental data, where it be- 
comes reduced by as much as a factor of 2 when passing 
from the 020 to 030 Brillouin zone. It also follows that the 
"waterfall" anomaly in the low-energy phonon dispersion 
is not specific to relaxor ferroelectrics, but it can be en- 
countered in other systems with interacting acoustic and 
soft optic branches, such as the ferroelectric perovskite 
BaTi0 3 . 



on the dispersion curves of several relaxor perovskites can 
be explained in the framework of coupled damped har- 
monic oscillator models describing the interaction of the 
acoustic and optic branches. This kink only appears on 
an apparent dispersion curve, obtained from maxima of 
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